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Abstract 

Let E = F(v) be the ground-state eigenvalue of the Schrodinger Hamil- 
tonian H = —A + vf(x), where the potential shape f(x) is symmetric and 
monotone increasing for x > 0, and the coupling parameter v is positive. If 
the kinetic potential f(s) associated with f(x) is denned by the transformation 
f(s) = F'(v), s = F(v) — vF'{y), then / can be reconstructed from F by the 
sequence /[ n+1 l = f o o /N. Convergence is proved for special classes of po- 

tential shape; for other test cases it is demonstrated numerically. The seed potential 
shape /'°] need not be 'close' to the limit /. 
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1. Introduction 

We consider the Schrddinger operator 

H = -A + vf(x) (1.1) 

defined on some suitable domain in L 2 (IR). The potential has two aspects: an 
'attractive' potential shape f(x), and a coupling parameter v > 0. We assume that 
f(x) is symmetric, non-constant, and monotone increasing for x > 0. Elementary 
trial functions can be designed for such an operator to prove that for each v > 
there is a discrete eigenvalue E = F(v) at the bottom of the spectrum. If, in 
addition to the the above properties, we also assume that / is continuous at x = 
and that it is piecewise analytic, then we are able to prove [1] that / is uniquely 
determined by F. The subject of the present paper is the reconstruction of the 
potential shape / from knowledge of the 'energy trajectory' F. This is an example 
of what we call 'geometric spectral inversion' [1, 2]. 

Geometric spectral inversion should be distinguished from the 'inverse problem 
in the coupling constant' which has been analysed in detail by Chadan et al [3-8]. 
In the latter problem the discrete part of the input data consists of the set {vi} of 
values of the coupling that yield a given fixed energy E. Inversion from excited- 
state energy trajectories Fk(v), k > 0, has also been studied, by a complete 
inversion of the WKB approximation for bound states [9]. For the ground-state 
energy trajectory F(y) = Fq(v) a constructive numerical inversion algorithm has 
been devised [10], and an inversion inequality has been established [11]. The work 
reported in the present paper also concerns inversion from the ground-state energy 
trajectory, but the approach uses functional methods which have a natural extension 
to the higher energy trajectories. 

Geometry is involved with this problem because we deal with a family of oper- 
ators depending on a continuous parameter v. This immediately leads to a family 
of spectral manifolds, and, more particularly, to the consideration of smooth trans- 
formations of potentials, and to the transformations which they in turn induce on 
the spectral manifolds. This is the environment in which we are able to construct 
the following functonal inversion sequence that is the central theme of the present 
paper: 
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f[n+l] = j Q f[n] 1 Q f[n] _ j Q R [n] _ 

A kinetic potential is the constrained mean value of the potential shape f(s) = </>, 
where the corresponding mean kinetic energy s = < — A> is held constant. It turns 
out that kinetic potentials may be obtained from the corresponding energy trajec- 
tory F by what is essentially a Legendre transformation [12] / <-> F given [13] 
by 



{f(s) = F'(v), s = F(v)—vF'(v)} <-> {F(v)/v = f(s)-sf'(s), l/v=-f'(s)}. 

(1.3) 

As we shall explain in more detail in Section 2, these transformations are well defined 
because of the definite convexities of F and /; they complete the definition of the 
inversion sequence (1.2), up to the choice of a starting seed potential f^(x). They 
differ from Legendre transformations only because of our choice of signs. The choice 
has been made so that the eigenvalue can be written (exactly) in the semi-classical 
forms 



E = F(y) = mm{s + vf(s)} = mm ^K^(x) + (1.4) 

where the kinetic- and potential-energy terms have the 'usual' signs. 

After more than 70 years of QM (and even more of the Sturm-Liouville prob- 
lem) it may appear to be an extravagance to seek to rewrite the min-max charac- 
terization of the spectrum in slightly different forms, with kinetic potentials and 
K functions. The main reason for our doing this is that the new representations 
allows us to tackle the following problem: if g is a smooth transformation, and we 
know the spectrum of —A + vp°\x), what is the spectrum of —A + vg(p°\x)) ? 
In the forward direction (obtaining eigenvalues corresponding to a given potential), 
an approximation called the 'envelope method' has been developed [13, 14]. The 
inversion sequence (1.2) was arrived at by an inversion of envelope theory, yielding 
a sequence of approximations for an initially unknown transformation g satisfying 
f(x) = g(fW(x)). 

In order to make this paper essentially self contained, the representation ap- 
paratus is outlined in Section 2. In Section 3 we use envelope theory to generate 
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the inversion sequence. In Section 4 it is proved that the energy trajectory for a 
pure power potential is inverted from an arbitrary pure-power seed in only two 
steps: thus = f in these cases. In Section 5 we consider the exactly soluble 
problem of the sech-squared potential f(x) = — sech 2 (x). Starting from the seed 
f^(x) = — 1 + x 2 /20, we are able to construct the first iteration /W exactly; we 
then continue the sequence by using numerical methods. This illustration is inter- 
esting because the seed potential shape is very different from that of the target 
f(x) and has a completely different, entirely discrete, spectrum. We consider also 
another sequence in which the seed is f^(x) = — 1/(1 + x/5). Convergence, which, 
of course, cannot be proved with the aid of a computer, is strongly indicated by 
both of these examples. 



2. Kinetic potentials and K functions 

The term 'kinetic potential' is short for 'minimum mean iso-kinetic potential'. 
If the Hamiltonian is H = — A+vf(x), where f(x) is potential shape, and T>(H) C 
L 2 (IR) is the domain of if, then the ground-state kinetic potential f(s) = fo(s) 
is defined [13, 14] by the expression 



f(s)= inf ty,/V0. (2-1) 

(V>,-AV0=s 

The extension of this definition to the higher discrete eigenvalues (for v sufficiently 
large) is straightforward [13] but not explicitely needed in the present paper. The 
idea is that the min-max computation of the discrete eigenvalues is carried out in 
two stages: in the first stage (2.1) the mean potential shape is found for each fixed 
value of the mean kinetic energy s; in the second and final stage we minimize over 
s. Thus we have arrive at the semi-classical expression which is the first equality 
of Eq.(1.4). It is well known that F(v) is concave (F"(v) < 0) and it follows 
immediately that f(s) is convex. More particularly, we have [2] 

F"{v)f" (a) = ~. (2.2) 

Thus, although kinetic potentials are defined by (2.1), the transformations (1.3) 
may be used in practice to go back and forth between F and /. 
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Kinetic potentials have been used to study smooth transformations of potentials 
and also linear combinations. The present work is an application of the first kind. 
Our goal is to devise a method of searching for a transformation g, which would 
convert the initial seed potential f^(x) into the (unknown) goal f(x) = g(f^). 
We shall summarize briefly how one proceeds in the forward direction, to approxi- 
mate F, if we know f(x). The K functions are then introduced, by a change of 
variable, so that the potential f(x) is exposed and can be extracted in a sequential 
inversion process. 

In the forward direction we assume that the lowest eigenvalue F^(v) of = 
—A + vf^(x) is known for all v > and we assume that f(x) is given; hence, 
since the potentials are symmetric and monotone for x > 0, we have defined the 
transformation function g. 'Tangential potentials' to g(f^) have the form a + 
&/ 0*0) where the coefficients a(t) and b(t) depend on the point of contact x = t 
of the tangential potential to the graph of f(x). Each one of these tangential 
potentials generates an energy trajectory of the form F(v) = av + F^(bv), and 
the envelope of this family (with respect to t) forms an approximation F A (v) to 
F(v). If the transformation g has definite convexity, then F A {v) will be either an 
upper or lower bound to F(v). It turns out [13] that all the calculations implied by 
this envelope approximation can be summarized nicely by kinetic potentials. Thus 
the whole procedure just described corresponds exactly to the expression: 

f"f A = gop \ (2.3) 

with m being replaced by an inequality in case g has definite convexity. Once 
we have an approximation f A , we immediately recover the corresponding energy 
trajectory F A from the general minimization formula (1.4). 

The formulation that reveals the potential shape is obtained when we use x 
instead of s as the minimization parameter. We achieve this by the following general 
definition of x and of the K function associated with / : 

f(x) = f(s), K^(x) = f-\f(x)). (2.4) 

The monotonicity of f{x) and of / guarantee that x and K are well defined. 
Since f~ 1 (f) is a convex function of /, the second equality in (1.4) immediately 
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follows [14]. In terms of K the envelope approximation (2.3) becomes simply 

tfW«tf[/ [01 ]. (2.5) 

Thus the envelope approximation involves the use of an approximate K function 
that no longer depends on /, and there is now the possibility that we can invert 
(1.4) to extract an approximation for the potential shape. 

We end this summary by listing some specific results that we shall need. First 
of all, the kinetic potentials and K functions obey [13, 14] the following elementary 
shift and scaling laws: 



f{x) - a + bf(x/t) =► |/» -> a + bf(st 2 ), K^{x) -> ItfW (|) } . (2.6) 
Pure power potentials are important examples which have the following formulas: 



f(x) = \x\« 




K(x) = - 



(2.7) 



where, if the bottom of the spectrum of —A + \x\ q is E(q), then the P numbers 
are given [14] by the folowing expressions with n = : 



(2+9) 

P n {q) = \E n {q)\ — 



2 + q 



2 + q 



q^O. 



(2.8) 



We have allowed for q < and for higher eigenvalues since the formulas are es- 
sentially the same. The P n (q) as functions of q are interesting in themselves [14]: 
they have been proved to be monotone increasing, they are probably concave, and 
P n (0) corresponds exactly to the log potential. By contrast the E n (q) are not so 
smooth: for example, they have infinite slopes at q = 0. But this is another story. 
An important observation is that the K functions for the pure powers are all of 
the form (P(q)/x) 2 and they are invariant with respect to both potential shifts 
and multipliers: thus a + b\x\ q has the same K function as does \x\ q . For the 

harmonic oscillator P n (2) = (n + i) 2 , n = 0,1,2, Other specific examples 

may be found in the references cited. 
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The last formulas we shall need are those for the ground state of the sech- 
squared potential: 

f(x) = -sech 2 (x) (/» = ^ , K(x) = sinh- 2 (2x)l . (2.9) 

I (s + S 1 ) 2 + s J 

3. The inversion sequence 

The inversion sequence (1.2) is based on the following idea. The goal is to find 
a transformation g so that / = g o p K We choose a seed /M, but, of course, / 
is unknown. In so far as the envelope approximation with as a basis is 'good', 
then an approximation f or g would be given by / = g^ o p°\ Thus we have 

g*gW=fofM-\ (3.!) 
Applying this approximate transformation to the seed we find: 

/ » fW = g W o /M = / o /M" 1 o /M = /o^l. (3.2) 

We now use /t 1 ! as the basis for another envelope approximation, and, by repetition, 
we have the ansatz (1.2), that is to say 

f[n+i] = fofWr 1 o /N = /oifW. (3.3) 

A useful practical device is to invert the second expression for F given in (1.4) 
to obtain 

K [f] (x) = max{F(v) -vf(x)}. (3.4) 

v>0 

The concavity of F(v) explains the max in this inversion, which, as it stands, is 
exact. In a situation where / is unknown, we have / on both sides and nothing 
can be done with this formal result. However, in the inversion sequence which we 
are considering, (3.4) is extremely useful. If we re-write (3.4) for stage [n] of the 
inversion sequence it becomes: 
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KM ( x ) = max (fW (v) - vf^ (x) ) . (3.5) 

v>0 I J 

In this application, the current potential shape p n * and consequently F^iy) can 
be found (by shooting methods) for each value of v. The minimization can then be 
performed even without differentiation (for example, by using a Fibonacci search) 
and this is a much more effective method for = /t n ' o /t n l than finding 

/t n '(s), finding the functional inverse, and applying the result to p n K 

4. Inversion for pure powers 

We now treat the case of pure-power potentials given by 

f(x) = A + B\x\ q , q>0, (4.1) 

where A and B > are arbitrary and fixed. We shall prove that, starting from 
another pure power as a seed, the inversion sequence converges in just two steps. 
The exact energy trajectory F(v) for the potential (4.1) is assumed known. Hence, 
so is the exact kinetic potential given by (2.7) and the general scaling rule (2.6), 
that is to say 

f(s) = A + B(^pj q . (4.2) 
We now suppose that a pure power is also used as a seed, thus we have 

fW(x) = a + b\x\ p KW(x)= (^^j , (4-3) 

where the parameters a, b > 0, p > are arbitrary and fixed. The first step of 
the inversion (1.4) therefore yields 

/W W = (/»A-B) (x) = A + B (^f 1 )". (4.4) 

The approximate potential f^(x) now has the correct x power dependence but 
has the wrong multiplying factor. Because of the invariance of the K functions to 
multipliers, this error is completely corrected at the next step, yielding: 
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K[11 ^ = (?) ^ f [2 \x) = [foK^{x)=A + B\x\ q . [4.5] 

This establishes our claim that power potentials are inverted without error in exactly 
two steps. 

The implications of this result are a little wider than one might first suspect. 
If the potential that is being reconstructed has the asymptotic form of a pure 
power for small or large x, say, then we know that the inversion sequence will very 
quickly produce an accurate approximation for that part of the potential shape. 
More generally, since the first step of the inversion process involves the construction 
of K^°\ the general invariance property K^ a+b ^ = K^* given in (2.6) means that 
the seed potential may be chosen without special consideration to gross features 
of / already arrived at by other methods. For example, the area (if the potential 
has area), or the starting value f(0) need not be incorporated in p°\ say, by 
adjusting a and b. 



5. More general examples 

We consider the problem of reconstructing the sech-squared potential f(x) = 
— sech 2 (x). We assume that the corresponding exact energy trajectory F(v) and, 
consequently, the kinetic potential f(s) are known. Thus [14]: 

/(.r) = cl.-'(.r) > { F(v) = -((v+l)i-l) , /» = 2 ' s 



4' 2J (s + s 2 )^ +s 

(5.1) 

We study two seeds. The first seed is essentially x 2 , but we use a scaled version of 
this for the purpose of illustration in Fig.(l). Thus we have 

This potential generates the exact eigenvalue 



(5.2) 
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which, like the potential itself, is very different from that of the goal. After the first 
iteration we obtain 



A graph of this potential is shown as fl in Fig.(l). In order to continue analytically 



be done numerically, with the aid of the inversion formula (3.5) for K. The first 5 
iterations shown in Fig.(l) suggest convergence of the series. 

As a second example we consider the initial potential given by 

/|0 ' W = -T+W5- (5 ' 4) 

In this case none of the steps can be carried out exactly. In Fig. (2) we show the first 
five iterations. Again, convergence is indicated by this sequence, with considerable 
progress being made in the first step. 

The numerical technique needed to solve these problems is not the main point 
of the present work. However, a few remarks about this are perhaps appropriate. As 
we showed in Ref.[10], by looking at three large values of v we can determine the 
best model of the form f(x) = A + B\x\ q that fits the given F(v) for small x < x a . 
We have used this method here with v = 10000 x {1, |, |} and x a = 0.2 in all cases. 
As indicated above, the inversion (3.5) was used to determine each function 
from the corresponding F^ n \ For all the graphs shown in the two figures, the range 
of v still to be considered for x > x a turned out to be 0.0008 < v < 175. With 
40 points beyond x a on each polygonal approximation for f^(x), a complete 
set of graphs for one figure took about 4| minutes to compute with a program 
written in C++ on a PentiumPro running at 200MHz. The exact expression for 
/M arising from the harmonic-oscillator starting point was very useful for verifying 
the behaviour of the program. 




(5.3) 



we would need to solve the problem with Hamiltonian = — A+vf^(x) exactly 
to find an expression for F^(v). We know no way of doing this. However, it can 
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6. Conclusion 

Progress has certainly been made with geometric spectral inversion. The results 
reported in this paper suggest strongly that in some suitable topology, the inversion 
sequence (1.4) converges to /. The 'natural' extension of (1.4) to excited states 
leads to the conjecture that each of the following sequences converges to / : 

f^/^f'o/^/^f 1 . k = 0,1,2,... (6.1) 

For the examples studied by the inversion of the WKB approximation, inversion 
improved rapidly as k increased and the view of the problem became more 'classi- 
cal'. 

If an energy trajectory F(v) is derived from a potential which vanishes at 
large distances, is bounded, and has area, it is straightforward [1] to 'normalize' 
the function F(v) by scaling so that it corresponds to a potential with area 2 and 
lowest point —1. The graphs of F (v) for normalized potentials with square-well, 
exponential, and sech-squared shapes look very similar: for small v they are asymp- 
totically like — v 2 , and for large v they satisfy lim„^ 00 {Fo(i')/'i>} = —1. These 
asymptotic features are the same for all such normalized potentials. We now know 
that encoded in the details of these F (v) curves for intermediate values of v are 
complete recipes for the corresponding potentials. If, as the WKB studies strongly 
suggest, the code could be unravelled for the excited states too, the situation would 
become very interesting. What this would mean is that, given any one energy tra- 
jectory Fkiy), we could reconstruct from it the underlying potential shape f(x) 
and then, by solving the problem in the forward direction, go on to find all the 
other trajectories {Fj(v)}j^k, and all the scattering data. For large k this would 
imply that a classical view alone could determine the potential and hence all the 
quantum phenomena which it might generate. 
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Figure (1) The energy trajectory F for the sech-squared potential f(x) = 
— sech 2 (x) is approximately inverted starting from the seed f^(x) = — 1 + x 2 /20. 
The first step can be completed analytically yielding fl = f^(x) = — 2/{l + 
a/1 + Ax 2 }. Four more steps {fk = /'^}| =2 of the inversion sequence approaching 
/ are performed numerically. 
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Figure (2) The energy trajectory F for the sech-squared potential f(x) = 
— sech 2 (x) is approximately inverted starting from the seed fO = f^(x) = —1/(1 + 
x/5). The first 5 steps {fk = f^}\ =1 of the inversion sequence approaching / 
are performed numerically. 



